# graficar histograma por area de poligonos juntos y por cobertura
fig,(ax1, ax2) = plt.subplots(nrows=1, ncols=2, figsize=(10, 7))
ax1.hist(
gdf_cob.Area_ha,
bins='sturges' # numero de agrupaciones se puede poner un numero o un metodo: 'auto', 'fd', 'doane', 'scott', 'stone', 'rice', 'sturges', or 'sqrt'
)
# agregar linea punteada con valor medio de area de poligonos
ax1.axvline(
x=gdf_cob.Area_ha.mean(),
color = 'red', #color de linea
linestyle='dashed', #estilo de linea
linewidth= 2, # ancho de linea
label='media: {:5.2f}'.format(gdf_cob.Area_ha.mean()) # label para linea punteada
)
ax1.set_xlabel('Área (ha)') # nombre eje x
ax1.set_ylabel('Frecuencia') # nombre eje y
ax1.set_title('Histograma de frecuencias')
ax1.legend(loc=0)
# graficar histograma por area de poligonos y categoria natural o antropica
ax2.hist(
[gdf_cob[gdf_cob['categ_cob']=='Antropica'].Area_ha, gdf_cob[gdf_cob['categ_cob']=='Natural'].Area_ha], # datos para histograma, van entre corchetes para graficar dos variables
bins='sturges', # numero de agrupaciones se puede poner un numero o un metodo: 'auto', 'fd', 'doane', 'scott', 'stone', 'rice', 'sturges', or 'sqrt'
label=['Antropico', 'Natural'] # label para barras de histograma, van en el mismo orden de las variables a graficar
)
# agregar linea con valor medio de area de poligonos natural
ax2.axvline(
x=gdf_cob[gdf_cob['categ_cob']=='Natural'].Area_ha.mean(), # valor de linea punteada
color = 'green', # color de linea
linestyle='dashed', # estilo de linea
linewidth= 1, # ancho de linea
label='media naturales: {:5.2f}'.format(gdf_cob[gdf_cob['categ_cob']=='Natural'].Area_ha.mean()) # label para leyenda
)
# agregar linea con valor medio de area de poligonos con actividad antropica
ax2.axvline(
x=gdf_cob[gdf_cob['categ_cob']=='Antropica'].Area_ha.mean(), # valor de linea punteada
color = 'blue', # color de linea
linestyle='dashed', # estilo de linea
linewidth= 1, # ancho de linea
label='media antropicos: {:5.2f}'.format(gdf_cob[gdf_cob['categ_cob']=='Antropica'].Area_ha.mean()) # label para leyenda
)
ax2.set_xlabel('Área (ha)') # nombre del eje x
ax2.set_ylabel('Frecuencia') # nombre del eje y
ax2.set_title('Histograma por categoria Antropica y Natural')
ax2.legend() # activa la leyenda para que se muestre en la grafica
fig.tight_layout() # ajusta las graficas para que no se sobrepongan la de la izquida con la de la derecha
fig.show()
